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v  Abstract 

Unsteady,  thin  layer  Navier-Stokes  equations  written 
in  rotor  coordinates  are  solved  using  a  partially  flux-split, 
implicit  numerical  algorithm  to  calculate  the  flowfields  of  a 
hovering  rotor  blade  at  subsonic  and  transonic  conditions. 
Numerical  results  are  in  good  agreement  with  experimental 
data  for  both  nonlifting  and  lifting  rectangular  blades.  For 
the  lifting  case,  the  rotor  wake  effects  are  modeled  by  ap¬ 
plying  a  correction  to  the  geometric  angle  of  attack  of  the 
blades.  Alternate  methods  of  calculating  the  hovering  rotor 
flowfields  in  blade-fixed  mode  that  have  the  same  circula¬ 
tion  distribution  as  hovering  blade  are  explored.  All  of  the 
results  presented  in  this  paper  are  computed  on  a  CRAY2 
supercomputer. 

l\ 

'  Introduction 

The  need  to  accurately  calculate  the  flowfield  of  a  he¬ 
licopter  rotor  in  both  hover  and  forward  flight  is  of  great 
practical  importance.  Unlike  the  flowfield  of  a  fixed  wing, 
the  flowfield  of  a  helicopter  rotor  is  generally  more  complex 
to  analyze  because  it  provides  some  of  the  most  complex 
challenges  to  be  found  in  the  field  of  applied  aerodynamics. 
This  complexity  stems  from  several  peculiar  problems  that 
are  unique  to  the  helicopter  rotor,  viz.,  a  radially  increasing 
blade  speed  that  is  responsible  for  a  high  concentration  of 
bound  circulation  over  the  outer  portion  of  the  blade  re¬ 
sulting  in  a  strong  trailed  vortex,  a  spiralling  wake  vortex 
sheet  remaining  initially  close  to  the  rotor  causing  strong 
blade-vortex  interactions,  a  high  centrifugal  force  field  in 
which  the  blades  operate,  a  relatively  large  steady  state 
out-of-plane  displacement  of  the  rotor  blades  and  aeroelas- 
tic  response  of  the  rotor  itself,  and  finally,  mutual  interac¬ 
tion  of  flowfields  of  main  rotor,  tail  rotor  and  the  fuselage. 
These  flowfields  are  often  characterized  by  transonic  con¬ 
ditions  and  associated  shock  waves  which  makes  the  flow 
more  susceptible  to  three  dimensionality  and  unsteadiness. 

The  operating  characteristics  of  such  rotory  wing  ve¬ 
hicles  are  strongly  influenced  by  the  the  vortex  wake.  The 
interaction  of  this  wake  with  the  following  blades  is  a  po¬ 
tential  source  of  noise  and  vibration  at  low  and  moderate 
flight  speeds.  Accurate  prediction  of  the  vortical  wake  is 
probably  the  most  important,  most  studied  and  the  most 
difficult  aspect  of  helicopter  flowfield.  Current  methods  of 
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analysis  of  the  wake  range  in  complexity  from  relatively 
simple  momentum-theory  applications  to  free  wake  lifting 
surface  methods.  In  between  these  extremes,  there  are  a  va¬ 
riety  of  so-called  prescribed-wake  models,  which  generally 
rely  on  some  degree  of  empiricism  to  determine  the  posi¬ 
tion  of  the  wake  vortices;  then  the  Biot-Savart  law  is  used 
to  calculate  the  induced  velocity  field  due  to  these  vortices 
at  the  plane  of  the  rotor  blades.  Although  such  models 
are  widely  used  in  current  prediction  techniques,  they  suf¬ 
fer  from  the  limitation  that  the  emipirical  determination 
of  the  wake  shape  ignores  some  of  the  important  details  of 
the  flowfield  such  as  the  mutual  interaction  between  various 
vortex  elements.  Further,  they  are  unreliable  for  unusual 
blade  pianforms  and/or  twist  distributions  which  are  often 
the  case  with  modern  helicopter  blade  shapes. 

The  current  thrust  in  calculating  the  rotor  flowfield 
more  accurately,  including  the  wake  effects,  is  pushing  the 
use  of  state-of-art  computational  fluid  dynamics  (CFD) 
codes  to  shed  more  light  on  the  understanding  of  this  prob¬ 
lem.  Finite  difference  codes  for  nonlinear  compressible  po¬ 
tential  equations  (Refs.  1-6)  and  the  Euler  equations  (Refs. 
7-14)  have  been  used  to  calculate  the  rotor  flowfields.  ini¬ 
tially  <!•  ■  i  loped  methods  using  the  potential  flow  and  the 
Euler  lormulations  were  primarily  limited  to  calculating 
nonliiting  rotor  flows  because  of  the  inherent  limitation  of 
not  being  able  to  model  the  vortex  wake  with  these  equa¬ 
tions,  although  the  Euler  formulation  has  in  it  the  neces¬ 
sary  physics  to  model  vorticity  transport  correctly.  These 
equations  basically  lack  the  physical  mechanism  needed  to 
generate  the  vortex  wake.  However,  in  conjunction  with 
wake  models,  such  as  CAMRAD  (Ref.  15)  and  HOVER 
(Ref.  16),  both  potential  flow  and  Euler  codes  have  been 
used  extensively  to  calculate  the  lifting  rotor  flowfields. 
The  standard  experimental  data  that  is  used  in  validat¬ 
ing  most  of  these  codes  has  been  the  two-bladed  rotor  data 
of  Caradonna  and  Tung  (Ref.  17).  An  excellent  review 
of  some  of  the  currently  available  inviscid  finite-difference 
numerical  methods  has  been  presented  by  Caradonna  and 
Tung  in  Ref.  18. 

As  mentioned  above,  tip  vortices  are  an  important  part 
of  the  helicopter  rotor  flow  field.  These  vortices,  which  are 
generated  at  the  tips  of  the  rotating  blades,  along  with 
the  helical  wake  vortex  sheet  have  tremendous  influence 
on  the  operating  characteristics  of  the  rotor.  Some  of  the 
common  practical  problems  caused  by  such  concentrated 
trailing  vortices  are  the  rotor  vib  ation  due  to  unsteady 


lift  fluctuation,  increased  induced  drag  and  the  an¬ 
noying  ’blade-slap’,  an  impulsive  noise  characteristic  are 
a  few  to  name.  Many  studies  have  been  made  to  date 
to  understand  and  reduce  the  influence  of  tip  vortices  by 
means  of  modifications  to  the  tip  geometry  of  the  rotating 
blades.  Various  analytical  and  numerical  studies  conducted 
are  basically  inviscid  in  nature  and  therefore  preclude  the 
mechanics  of  the  physics  necessary  to  model  correctly  the 
formation  of  the  tip  vortex  which  involves  the  the  com¬ 
plex  three-dimensional  viscous  flowfield  in  the  tip  region. 
The  thin  layer  Navier-Stokes  simulations  of  tip-flows  have 
been  attempted  only  recently  after  bigger  and  faster  super¬ 
computers  became  available.  The  studies  of  Mansour  (Ref. 
19)  and  Srinivasan  et  al  (Ref.  20)  show  limited  success  in 
simulating  the  complex  tip  flows.  Simulation  of  the  com¬ 
plete  vortex  wake  now  appears  possible  with  proper  Navier- 
Stokes  algorithm  in  conjunction  with  patched  and/or  zonal 
grid  topology  to  discretize  the  flowfield.  Understanding  the 
mechanism  of  the  formation  process  of  the  tip  vortex  and 
its  subsequent  roll-up  would  provide  a  proper  insight  to 
modify  these  tip  flows  and  alleviate  some  of  the  problems 
caused  by  them.  The  ability  to  preserve  the  concentrated 
vortices  in  the  finite  difference  grid  without  numerical  dif¬ 
fusion  (Ref.  21)  has  been  the  biggest  set  back  until  now 
for  much  progress  in  this  area.  Even  the  most  advanced 
computational  techniques,  that  use  spatial  central  differ¬ 
encing,  lack  proper  mechanism  to  preserve  concentrated 
tip  vortices  and  convect  them  in  the  flowfield  without  nu¬ 
merical  diffusion.  However,  the  recently  developed  upwind 
schemes  in  conjunction  with  a  proper  grid  choice  appear 
very  promising  to  preserve  and  convect  concentrated  vor¬ 
tices.  Alternatively,  if  the  properly  captured  tip  vortex  is 
analytically  represented,  then  prescribed  vortex  methods 
(Refs.  22-23)  could  be  applied  to  calculate  the  vortex  wake 
development  for  several  rotations  of  the  blade.  These  meth¬ 
ods  have  demonstrated  the  ability  to  preserve  and  convect 


fixed  blade),  keeping  the  circulation  distribution  the  same 
as  that  of  a  hovering  blade,  are  explored  and  the  results 
for  these  are  compared  with  the  hover  results.  Governing 
equations  and  numerical  method  are  presented  in  the  next 
section  followed  by  results  and  concluding  remarks. 

Governing  Equations  and  Numerical  Scheme 

The  governing  partial  differential  equations  are  the 
unsteady, thin-layer  Navier-Stokes  equations.  For  general¬ 
ity,  these  equations  are  transformed  from  the  Cartesian 
reference  frame  to  the  arbitrary  curvilinear  space  ({,  q,  f, 
r)  while  retaining  strong  conservation  law-form  to  capture 
shock  waves.  The  transformed  equations  are  given  by  (Ref. 
25) 


drQ  +  d^E  +  dvF  +  d<G  =  Re"‘dfS  (1) 
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concentrated  vortices  even  in  very  coarse  grid  regions  with¬ 
out  significant  numerical  diffusion. 

The  use  of  Navier-Stokes  codes  to  model  the  rotor  flow- 
fields  have  been  limited  in  the  past  primarily  because  of  the 
large  computer  memory  and  CPU  time  requirements.  In 
fact,  these  codes  may  not  be  very  much  more  expensive  to 
run  than  some  of  the  Euler  codes,  but  to  have  a  meaning¬ 
ful  flow  definition  in  the  tip  region  and  in  the  wake  might 
make  such  a  procedure  very  expensive.  Recently,  Wake  and 
Sankar  (Ref.  24)  have  presented  some  nonlifting  and  lifting 
calculations  for  a  rotor  in  hover  using  a  coarse  grid  with 
a  poor  definition  of  the  tip  region.  The  results  for  surface 
pressure  are  in  fair  agreement  with  experimental  data.  The 
lifting  calculations  used  a  correction  to  the  geometric  an¬ 
gle  of  attack  of  the  blade  obtained  from  the  lifting  surface 
code  of  Ref.  16.  In  the  present  study,  using  a  good  defi¬ 
nition  of  the  tip  region,  thin-layer  Navier-Stokes  equations 
are  solved  for  the  flowfield  of  a  rectangular  blade  in  hover 
with  a  view  to  capture  the  tip  vortex  also.  Both  nonlift¬ 
ing  and  lifting  cases  have  been  calculated  with  subcritical 
and  supercritical  tip  Mach  numbers.  Surface  pressure  dis¬ 
tributions  and  tip  flow  data  are  presented  and  compared 
with  experimental  data.  Alternate  methods  of  calculating 
the  hover  flowfield  in  a  non-hover  mode  (like  an  isolated 


The  viscous  flux  vector  5,  written  here  in  the  limit  of  thin- 
layer  approximation,  is  given  by 
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The  primitive  variables  of  Eq,  (1),  viz.,  the  density 
p,  the  mass  fluxes  pu,pv,pw  and  the  energy  per  unit  volume 
e,  are  normalized  by  the  free-stream  reference  quantities. 
The  reference  length  and  velocity  scales  are  the  chord  of 
the  rotor  blade  and  free  stream  speed  of  sound  respectively. 
Other  nondimensional  quantities  appearing  in  the  above 
equations  are  the  Reynolds  number  Rt  and  the  Prandtl 


•number  Pr\  p  is  the  dynamic  viscosity.  The  relations  for 
the  contravariant  velocities  U,  V  and  W ,  the  Jacobian  of 
transformation  J,  and  the  metrics  of  the  transformation 
( 6 M .t?. )  *"d  (ft ,  f, , ft,, f, )  can  be  found 
in  Ref.  2S. 

The  velocity  components  u,v,w  and  the  pressure, 
p,  are  related  to  the  total  energy  per  unit  volume,  e,  through 
the  equation  of  state  for  a  perfect  gas  by 


p-h- 1)(*  -  ^(“3  +  1’2 


w*)) 


(5) 


where  *y  is  the  ratio  of  specific  heats. 

In  the  above  equations  u,v,  and  w  are  the  Carte¬ 
sian  components  of  the  velocity  in  the  intertial  coordinate 
system  (x,y, z,t).  In  the  present  formulation  Eq.  (1)  is 
solved  in  the  inertial  frame  of  reference.  The  intertial  co¬ 
ordinates  X  =  (x,y,z,t)  are  related  to  the  blade  fixed  co¬ 
ordinates  Xb  =  (i,  y,  £,t)  through  the  relation  given  by 


X(x,y,z)  =  R(t)Xb(i,y,i) 

t  -  t 

where  R(t)  is  the  rotational  matrix  (Ref.  26)  given  by 
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Here  n  is  the  reduced  frequency  of  the  rotor  and  Qt 
represents  the  azimuth  sweep  of  the  rotor  blade.  In  view  of 
this  relation,  Eq.  7,  the  metrics  in  inertial  reference  frame 
are  related  to  those  in  the  blade-fixed  frame  of  reference 
through 

=  (iCosHf  -  (jsinflt 
fn  =  fisin  Ot  +  ( jcos  nt 
f,  =  fi 
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The  equations  set,  Eq.  1,  is  solved  using  an  im¬ 
plicit,  approximately-factored  numerical  scheme  that  uses 
spatial  central  differencii.’  in  the  r>  and  f  directions  and 
upwind  differencing  in  the  f  direction  developed  by  Ying 
et  al  (Ref.  27).  The  flux  vector  F  has  been  split  in  to  F* 
and  F~  according  to  its  eigenvalues.  Artificial  dissipation 
terms  (second-  and  fourth-order)  have  been  added  in  the 
central  differencing  directions  for  stability  reasons  (see  Ref. 
25).  The  factored  operators  can  be  solved  by  sweeping  in 


the  {  direction  and  inverting  tridiagonal  matrices  with  5  x 
5  blocks  for  the  other  two  directions.  Currently,  significant 
part  of  the  computational  time  is  taken  to  form  the  plus 
and  minus  Jacobian  matrices  for  the  flux  vector  F  with 
this  numerical  scheme.  The  numerical  code  is  vectorized 
for  the  CRAY2  supercomputer. 

A  body  conforming  finite-difference  grid  has  been 
used  for  the  rectangular  blade  having  a  rounded-tip  cap 
and  consists  of  warped  spherical  0-0  grid  topology.  The 
flowfield  grid  is  numerically  generated  using  the  three- 
dimensional  hyperbolic  grid  solver  of  Steger  and  Chaussee 
(Ref.  28)  with  proper  clustering  in  the  leading  and  trail¬ 
ing  edge  regions  and  in  the  tip  region.  The  grid  is  nearly 
orthogonal  at  the  surface  and  the  spacing  in  the  normal  di¬ 
rection  at  the  surface  is  chosen  to  be  0.00006  of  the  chord. 
All  the  computations  were  done  with  one  grid  topology  hav¬ 
ing  155  points  in  the  periodic  direction  around  the  airfoil, 
and  66  points  each  in  the  spanwise  and  normal  directions, 
for  a  total  of  about  700,000  grid  points.  The  grid  boundary 
is  chosen  to  be  at  10  chords  in  all  directions. 

The  boundary  conditions  consist  of  surface  bound¬ 
ary  conditions  and  farfield  boundary’  conditions  and  are  ap¬ 
plied  explicitly.  For  the  nonrotating  blade  the  noslip  con¬ 
dition  is  enforced  at  the  wall  by  setting  U,  V'  and  W'  to 
be  zero  and  £(,  It  and  ft  are  zero  as  the  grid  is  station¬ 
ary.  For  the  rotating  blade,  however,  U ,  V  and  W  are  still 
set  to  zero  but  £t,  rj,  and  ft  are  nonzero  as  the  blade  (and 
the  grid  attached  to  it)  is  moving  in  azimuth.  Also,  at  the 
wall  the  density  is  determined  by  assuming  adiabatic  wall 
condition.  The  pressure  along  the  body  surface  is  calcu¬ 
lated  from  the  normal  momentum  relation  (see  for  example 
Ref.  25).  Having  known  the  density  and  pressure,  the  total 
energy  is  determined  from  the  equation  of  state. 

At  the  farfield  boundary  the  flow  quantities  are  ei¬ 
ther  fixed  or  extrapolated  from  the  interior  depending  on 
whether  the  flow  is  subsonic  or  supersonic  and  if  it  is  of 
inflow-  or  outflow-type  at  the  boundary.  The  characteristic 
velocities  of  the  Euler  equations  determine  the  number  of 
flow  properties  to  be  specified  to  control  the  reflections  of 
waves  from  the  boundaries.  For  subsonic-inflcxv  boundary, 
four  quantities  must  be  specified.  Thus  density  is  extrapo¬ 
lated  while  the  velocities  and  the  total  energy  are  specified 
by  the  free  stream  values.  For  supersonic-inflow,  all  flow- 
quantities  are  specified.  At  subsonic-outflow  boundaries, 
only  one  quantity  is  specified,  viz.,  pressure  is  fixed.  For 
supersonic-outflow  condition  all  flow  quantities  are  extrap¬ 
olated  from  the  interior.  At  the  plane  containing  the  blade 
root  dQ/dy  =  0  is  imposed. 

Results  and  Discussion 

Both  time  accurate  and  time  asymptotic  (steady 
state)  calculations  are  performed  in  this  study  for  an  as¬ 
pect  ratio  6  rectangular  blade  having  no  twist  or  taper. 
The  blade,  which  has  a  round  tip-cap  (body  of  revolution), 
is  made  up  of  NACA0012  airfoil  section.  The  rotating  blade 
calculations  presented  here  have  been  chosen  to  correspond 
to  the  experimental  test  conditions  of  a  two-bladed  hov¬ 
ering  rotor  of  Caradonna  and  Tung  (Ref  17).  At  the 


Reynolds  number  corresponding  to  the  tip  speeds  in  this 
test,  the  boundary  layer  can  be  assumed  to  be  turbulent 
over  the  entire  blade  and  Baldwin  and  Lomax  algebraic 
turbulence  model  (Ref.  29)  is  used  to  calculate  the  tur¬ 
bulent  eddy  viscosity.  A  typical  solution,  with  vectorized 
code  for  CRAY2.  required  about  700-1000  marching  steps 
(approximate!)  45-60  degrees  of  azimuth  travel)  to  reach 
quasi-steady  flow  conditions  with  CPI'  time  per  time  step 
per  grid  point  of  8.5  x  1CT5  sec.  Time  accurate  calcula¬ 
tions  were  run  with  a  maximum  value  for  the  time  step  of 
At  =  0.01,  whereas  the  fixed  blade  option  used  varable  time 
step  option  (Ref.  21)  with  At  of  5  to  calculate  steady  state 
flowfield. 

Hovering  Blade 

As  mentioned  before,  the  rotating  blade  calcula¬ 
tions  presented  here  correspond  to  the  two-bladed  hover 
test  conditions  of  Caradonna  and  Tung  (Ref.  17).  Hover 
calculations  have  been  done  in  a  time  accurate  fashion.  In 
practice,  a  hovering  rotor  flowfield  is  quasi-steady  in  blade- 
fixed  coordinates  system.  Since  the  governing  equations  are 
being  solved  in  the  inertial  reference  frame  in  the  present 
approach,  the  flowfield  never  reaches  steady  state  in  this 
reference  frame  This  approach  was  preferred  since  it  is 
easily  extendable  to  the  forward  flight  case. 

Figure  1  shows  the  chordwise  surface  pressure  dis¬ 
tributions  at  several  radial  stations  for  the  nonlifting  rotor 
with  tip  Mach  number  A/(lf  -  0.52  and  a  corresponding 
Reynolds  number  of  2.32  million.  As  seen,  the  calculated 
results  are  in  excellent  agreement  with  experimental  data 
at  all  radial  stations.  Accurate  calculation  of  lifting  rotor 
flow-fields  is  possible  only  if  the  induced  effects  of  the  wake 
are  properly  included  in  the  analysis.  At  present  a  num¬ 
ber  of  methods  are  available  to  model  the  effects  of  the 
wake  of  a  hovering  rotor  All  these  methods  sec  .i  to  give 
approximately  the  same  kind  of  results  when  applied  to  a 
particular  problem,  not  dominated  by  viscous  effects,  as  is 
shown  in  the  review  of  Caradonna  and  Tung  (Ref.  18). 
With  this  observation,  a  uniform  correction  to  the  angle 
of  attack  of  the  blade  has  been  made  in  this  study  based 
on  the  estimates  of  induced  downwash  for  the  experimental 
test  rotor  configuration  given  by  Agarwal  and  Deese  (Ref. 
11)  Figure  2,  reproduced  from  Ref.  1 1 ,  shows  plots  of  sec¬ 
tional  induced  velocity  estimated  from  a  free- wake  analysis 
program  for  three  speeds  of  rotor  and  thrust  coefficients  of 
an  experimental  configuration  Over  a  section  of  the  blade, 
approximately  from  0  4  to  0  9  radius  in  this  figure,  the  in¬ 
duced  downwash  given  by  the  ratio  of  sectional  induced 
velocity  to  the  local  blade  speed  is  nearly  constant  equal  to 
-3  8  degrees  for  the  entire  range  of  test  conditions  In  the 
present  calculations  this  induced  downwash  is  chosen  as  a 
representative  value  and  is  assumed  constant  for  the  entire 
blade 

The  effective  pitch  of  the  hovering  blade  is  then  the 
difference  of  the  geometric  angle  of  attack  and  the  induced 
downwash  estimated  above  Using  this  estimate,  lifting  cal¬ 
culations  have  been  performed  for  hovering  blades  set  at  an 
effective  pitch  of  4  2  degrees  and  having  tip  Mach  numbers 
of  0  44  and  0  *>77.  respectively  1  ires  3  and  4  show  the 


computed  results  for  these  cases  in  the  form  of  surface  pres¬ 
sure  distributions  for  representative  blade  radial  stations 
compared  with  the  experimental  data  of  Caradonna  and 
Tung  (Ref.  17).  The  comparison  shows  very  good  agree¬ 
ment.  at  least  for  the  radial  stations  between  0.6  to  0.95 
of  radius,  for  both  subcriticai  and  supercritical  cases  The 
agreement  progressively  deteriorates  for  radial  stations  le«s 
than  0.6R  which  probably  s  expected  from  the  data  for 
Fig.  2.  Agarwal  and  Deese  (Ref.  11)  also  calculated  the 
same  flow  using  the  finite  volume  Euler  formulation  and  the 
same  induced  downwash  correction.  They  also  get  similar 
agreement  with  the  experimental  data  for  the  subcriticai 
case.  However,  the  shock  locations  are  over  estimated  for 
supercritical  case  in  their  calculations. 

Additional  flowfield  data  for  the  tip  region  are  pre¬ 
sented  in  Figs.  5-10  for  the  two  lifting  cases  calculated  here. 
It  is  well  known  that  the  formation  of  tip  vortex  involves 
complex  three-dimensional  flow  separation  in  the  tip  region 
as  a  necessary  condition.  Figure  5a  shows  a  view  of  the  sur¬ 
face  particle  flow  traces  in  the  tip  region  for  the  subcriticai 
lifting  configuration  corresponding  to  a  tip  Mach  number 
of  M i,r  =  0.44.  It  is  generated  by  releasing  fluid  particle 
tracers  at  one  grid  point  above  the  surface  and  confining 
the  flow  to  stay  in  that  plane.  This  view  is  supposed  to 
mimic  the  surface  oil  flow  pattern  often  used  in  laboratory 
experiments.  The  separation  and  reatlachment  lines  of  the 
flow  are  marked  by  the  symbols  S  and  R  respectively  in 
Fig.  5a.  The  extent  of  separation  on  the  upper  surface, 
inboard  of  tip.  is  much  larger  compared  to  the  mild  sep¬ 
aration  seen  on  the  lower  surface  of  the  tip  region  The 
braiding  of  flow  particle  tracers,  released  from  different  lo¬ 
cations  on  the  upper  and  lower  surfaces,  in  the  tip  region 
show  the  formation  process  of  tip  vortex  in  Fig.  6a.  In 
contrast  to  a  nonrotating  blade  (Ref.  20),  tne  braiding  of 
particle  tracers  from  upper  and  lower  surfaces  is  delayed 
until  after  the  vortex  lifts-off  from  the  surface  The  initial 
braiding  process  is  comprised  of  mostly  the  particles  from 
the  upper  surface  as  is  clearly  seen  in  Fig.  6a  The  strength 
of  the  vortex  is  determined,  as  before  (Ref.  20).  by  the  line 
integral  of  the  velocity  vector  over  a  closed  path  enclosing 
the  vortex.  Estimations  done  at  several  locations  in  the 
downstream  wake  by  this  method  gave  a  value  of  0.08  to 
C  09  depending  on  the  size  of  the  line  integral  path  chosen 
for  the  vortex  strength.  The  integrated  lift  from  the  blade 
pressure  distributions  was  found  to  be  0  19  The  vorticity 
contours  shown  in  Fig.  6a  show  the  size  and  shape  of  the 
tip  vortex  As  expected,  the  coarse  grid  has  smeared-off 
the  tightly  wound  vortex  as  seen  in  this  figure. 

Similar  results  are  presented  for  the  lifting  hov¬ 
ering  rotor  blade  with  A/,tp  -  0.877  in  Figs.  7-10.  The 
pressure  contours  of  Fig.  7a  show  the  extent  of  transonic 
.low  and  steeping  of  the  shock  wave  towards  the  tip  region. 
The  surface  particle  flow  pattern  of  Fig.  8a  clearly  iden¬ 
tify  the  regions  of  separation  and  reattachment.  As  before, 
the  separation  seen  on  the  lower  surface  in  the  tip  region 
is  milder  compared  to  the  extent  seen  inboard  of  the  up¬ 
per  surface.  The  formation  process  of  the  tip  vortex,  seen 
in  Fig  9a.  shows  the  braiding  of  particle  tracers  initially 
consisting  mostly  from  the  upper  surface  before  lifting-off 


the  upper  surface  inboard  of  the  tip.  Further  braiding  of 
this  from  the  particles  from  the  lower  surface  occurs  in  the 
wake  during  the  roll-up  process.  The  view  looking  at  the 
tip,  seen  in  Fig.  9a,  clearly  shows  the  lifting-off  of  the  tip 
vortex  from  the  separated  region  on  the  upper  surface.  A 
farfield  view  of  this  vortex  is  shown  in  Fig.  10a.  The  vor- 
tieity  contours  shown  in  this  figure  at  different  x-iocations 
in  the  wake  suggests  approximate  shape  of  the  vortex.  The 
strength  of  this  vortex  estimated,  as  before,  gate  a  value 
of  0.07-0.08.  The  integrated  value  of  the  lift  coefficient  for 
this  case  is  0.17. 


The  results  presented  in  the  above  section  for  the 
hovering  blade  were  computed  in  a  time  accurate  manner  as 
mentioned  before  and  these  are  time  consuming  and  expen¬ 
sive.  So  this  section  explores  alternate  methods  to  compute 
the  quasi-steady  flowfield  of  the  hovering  rotor,  in  particu¬ 
lar,  as  flowfield  of  a  fixed  blade  with  the  same  circulation 
distribution  as  that  of  the  hovering  rotor  blade  with  the 
same  tip  Mach  number.  Comparing  the  circulation  distri¬ 
bution  for  these  two  modes  (fixed  and  rotating  blade),  one 
can  immediately  come  up  with  two  different  ways  of  gener¬ 
ating  approximately  the  same  circulation  distribution  on  a 
fixed  blade  as  that  of  the  hovering  blade  keeping  the  plan- 
form  of  the  blade  same  by  a)  having  the  flow  Mach  number 
distributed  as  a  function  of  the  radial  distance  exactly  like 
a  hovering  blade,  keeping  every  thing  else  same,  and  6)  al¬ 
ternatively  keeping  the  flow  Mach  number  uniform  for  the 
entire  blade  equal  to  the  tip  speed  and  then  have  a  twist 
distribution  along  the  blade  which  decreases  from  the  tip 
to  the  root  of  the  blade.  This  means  that  the  fixed  blade 
will  have  variable  twist  increasing  towards  the  tip  to  a  value 
equal  to  the  value  of  effective  pitch  of  the  hovering  blade. 
A  third  option,  that  of  increasing  the  chord  linearly  from 
the  root  to  the  tip,  was  not  explored  in  this  investigation. 

With  the  above  reasoning,  steady  state  flowfieids 
were  calculated  just  like  a  fixed  isolated  blade  flowfield  with 
the  free  stream  conditions  as  discussed  above.  These  calcu¬ 
lations  used  variable  time  step  option  suggested  in  Ref.  21. 
Figures  11  and  12  show  these  results  in  the  form  of  surface 
pressure  distributions  compared  with  the  hovering  rotor  re¬ 
sults  for  both  subcritical  and  supercritical  cases.  The  sub- 
critical  results  presented  in  Fig.  11  show  surprisingly  very- 
good  agreement  with  the  hover  results  for  both  the  options 
of  variable  twist  and  variable  Mach  number.  Considering 
close  agreement  of  the  results  even  at  the  radial  station  near 
the  tip  of  the  blade,  the  influence  of  the  centrifugal  forces 
present  in  rotating  blade  seem  to  have  very  little  influence 
in  modifying  the  pressure  field  in  the  tip  region.  There 
are  very  small  differences  in  the  surface  particle  flow  pat¬ 
terns  and  the  locations  of  vortex  lift-off  from  the  surface  for 
these  cases  compared  to  the  hovering  blade  as  seen  in  Figs. 
6  and  7.  The  fixed-blade  configurations  seem  to  produce 
tightly  wound  vortex  even  before  leaving  the  blade,  but  its 
strength  appears  to  be  within  5-10  percent  of  the  hovering 
blade  value. 

The  supercritical  results,  presented  in  Fig.  12.  al¬ 
though  dominated  by  the  transonic  shocks  in  the  tip  region 


show  very  good  agreement  of  the  variable  Mach  number 
case  with  the  hover  case.  The  variable  twist  option  does 
not  seem  to  perform  as  well  in  the  transonic  regime.  This 
is  not  surprising  since  high  flow  Mach  number  (equal  to  lip 
Mach  number  of  rotor)  exists  all  along  the  span  for  this 
non-rotating  case.  Figures  7b  and  7c  show  the  upper  sur¬ 
face  pressure  contour  plots  for  these  cases  and  comparing 
these  with  that  of  hovering  blade  in  Fig  7a  demonstrates 
the  closeness  of  the  variable  Mach  number  fixed-blade  case 
with  the  hovering  case.  The  variable  twist  case  produces 
too  strong  a  shock  wave  along  the  entire  blade  with  the 
consequece  of  producing  a  totally  different  flowfield.  The 
close  agreement  of  the  flowfield  for  the  cases  of  Figs.  7a 
and  7b  suggests  that  even  at  this  supercritical  flow  condi¬ 
tion,  the  influence  of  the  centrifugal  forces  appear  to  have 
minimal  overall  effect  on  the  flowfield.  However,  the  sur¬ 
face  particle  flow  pattern  seen  in  Figs.  8a-8c  shows  differ¬ 
ences  in  the  flowfieids  in  the  tip  region  for  these  cases.  The 
rotating-blade  case  of  Fig.  8a  and  the  fixed-blade  case  of 
Fig. 8b  although  have  similar  surface  pressure  distributions, 
the  flow  in  the  tip  region  near  the  surface  appear  different 
The  difference  between  the  two  cases  should  come  from 
the  difference  in  here  force  fields,  which  means  that  the 
presence  of  centrifugal  force  in  the  rotating-blade  case  may- 
have  produced  a  small-scale  local  separation  on  the  lower 
surface  of  the  tip  in  addition  to  the  separation  seen  on  the 
upper  surface  of  the  blade.  This  separation  is  clearly  is  seen 
in  Fig.  8a.  The  flow  in  the  tip  region  for  the  fixed-blade 
cases  of  Figs.  8b  and  8c  are  largely  similar  even  though 
the  flowfield  for  the  case  of  Fig.  8c  is  dominated  by  the 
strong  shock  wave  and  its  associated  effects.  However,  the 
nearfield  views  of  the  tip  vortices  shown  in  Figs.  9a-9c  ap¬ 
pear  nearly  identical,  although  the  lift-off  of  the  tip  vortex 
for  hovering  blade  occurs  well  inboard  of  the  tip  on  the  up¬ 
per  surface  in  contrast  to  the  fixed  blade  cases  of  Figs.  9b 
and  9c  for  which  the  lift-off  appear  to  occur  right  in  the  tip 
region  on  the  upper  surface  of  the  blade.  The  farfield  views 
of  the  tip  vortices  for  these  cases  is  shown  in  Figs.  10a-  10c 
along  with  vorticity  contours  insert  showing  the  cross  sec¬ 
tions  of  the  tip  vortices  and  approximate  shapes  at  various 
distances  from  the  blade.  While  the  strength  of  the  tip  vor¬ 
tex  for  the  variable  Mach  number  case  is  almost  identical 
to  that  of  the  hovering  blade  case,  the  variable  twist  case 
produced  a  vortex  approximately  2.5  times  that  of  the  hov¬ 
ering  case.  The  vortex  shapes,  determined  by  the  vorticity 
contour  inserts,  appear  to  be  same  for  all  cases.  It  should 
be  emphasised  again  that  these  time  accurate  and  steady 
state  calculations  were  done  on  the  same  grid  topology  to 
remove  the  grid  dependency  from  the  comparisons. 

Conclusions 

Unsteady,  thin  layer  Navier-Stokes  equations  writ¬ 
ten  in  rotor  coordinates  are  solved  using  a  flux-split  approx¬ 
imately  factored,  implicit,  numerical  algorithm  to  calculate 
the  quasi-steady  flowfield  of  a  hovering  rotor  blade.  The 
test  cases  chosen  correspond  to  the  experimental  model 
hover  test  conditions  of  Caradonna  and  Tung  (Ref.  17). 
The  numerical  results  compare  very  well  with  the  experi¬ 
mental  data  for  both  nonlifting  and  lifting  cases.  The  in¬ 
duced  wake  effects  in  the  lifting  calculations  were  accounted 
as  a  correction  to  the  geometric  angle  of  attack  (pitch).  Al- 


lernate  methods  are  explored  to  calculate  the  hovering  ro¬ 
tor  flowfield  as  steady  state  flowfield  on  fixed  isolated  blade 
keeping  the  same  ciculation  distribution  as  that  of  the  hov¬ 
ering  blade.  Of  the  two  options  considered,  the  variable  free 
stream  Mach  number  case  gave  almost  identical  results  as 
that  of  a  rotor  at  both  subcritical  and  supercritical  flow 
conditions.  The  variable  twist  option,  on  the  other  hand, 
gave  similar  results  only  under  subcritical  flow  conditions; 
the  supercritical  flow  condition  was  dominated  by  strong 
transonic  shocks.  Under  conditions  where  the  fixed  blade 
flowfield  closely  agreed  with  that  of  hovering  blade,  the 
influence  of  the  centrifugal  forces  of  the  rotating  blade  ap¬ 
peared  to  have  minimum  influence  on  the  overall  flowfield 
properties.  While  these  conclusions  are  preliminary  and 
primarily  based  on  the  comparison  of  surface  pressures  and 
vortex  strength  estimates,  further  quantitative  comparison 
of  the  vortex  structure  is  needed  for  a  clear  understanding 
of  the  similarities  and  differences. 
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